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The Faddeev-Yakubowski equations have been solved in configuration space for the 
four nucleons system. Results for bound and scattering states in the isospin and S-wave 
approximation for different (T,S) channels are presented. The n-t elastic cross section 
has been also calculated with realistic interactions. A special interest was devoted to the 
description of the observed resonant structure at 7] a & « 3 MeV which is well reproduced 
only by the simple MT I-III potential. 

1. INTRODUCTION 

We present in this contribution the results obtained by solving the Faddeev-Yakubowski 
(FY) equations in configuration space for the 4N system with special interest in their scat- 
tering states. The A=4 chart presents a very rich continuum and constitutes a challenge 
for the theoretical models. We have limited our study to the energy region below the first 
inelastic channel. All the calculations have been carried out in the isospin approximation, 
i.e neglecting Coulomb forces. Next section will be devoted to a short presentation of 
the FY formalism. In Section 3 we will discuss the results: first those obtained in the 
S-wave approximation for the different (T,S) channels and secondly those concerning the 
n-t scattering with realistic interactions. 



2. FORMALISM 

The FY equations for 4 interacting particles are obtained by first splitting the total 
wavefunction \1/ in the usual Faddeev amplitudes, associated with each of the inter- 
acting pair. Each of these amplitudes is in its turn splitted in 3 parts, the FY amplitudes, 
corresponding to the different asymptotics of the remaining two non interacting particles. 
The Schrodinger equation is equivalent to the set of 18 coupled equations for which it is 
possible to define appropriate boundary conditions ensuring the unicity of the solution. 

In the case of identical particles, the 18 FY amplitudes can be obtained by the action of 
the permutation operators P^ on two of them, denoted K and H. The K — H amplitudes 
satisfy the following equations: 

(E-H -V)K = V[(P 23 + P n )(e + P u )K + e(P 23 + P 13 )H] (1) 
(E-H -V)H = V[(P 13 P 24 + P 14 P 23 )K + P 13 P 24 H] (2) 



in which e = ±1 depending on whether we are dealing with bosons or fermions. Each 
amplitude F — K,H is considered as function of its own set of Jacobi coordinates x, y, z 
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and expanded in angular momentum variables for each coordinate according to 

< xyz\F >= J2 Y a (x,y,z) (3) 

where Y a are generalized tripolar harmonics containing spin, isospin and angular mo- 
mentum variables and the functions F a , called FY components, are the unknowns. The 
label a holds for the set of intermediate quantum numbers defined in a given coupling 
scheme and includes the specification for the type of amplitudes K or H. We have used 
the following couplings: 



K amplitudes {[(^2)^3]^ U) 
H amplitudes (ii^)^^^ 




.hi 




Jxy ) J 

where t { and are the isospin and spin of the individual particles and (T, J) are respec- 
tively the total isospin and angular momentum of the four-body system. Each component 
F a is thus labelled by a set of 12 quantum numbers to which the symmetry proper- 
ties of the wavefunction impose the additional constraints: ^—\Y x+Tx+lx = £ for K and 
f_iy x +T X +ix = (_i\a-y+r y +i v _ £ f or g 

The boundary conditions are implemented by imposing at large enough values of z 
the asymptotic behaviour of the solutions. Indeed, at large values of z and for energies 
below the first inelastic threshold, the solution of equation ([D factorizes into a bound 
state solution of the 3N Faddeev equations and a plane wave in the z direction whereas 
the solution of (0) vanishes. The scattering observables are directly extracted from the 
logarithmic derivative of the K amplitude in the asymptotic region. The asymptotic 3N 
state is calculated with the same numerical scheme than the one used to solve the 4-body 
problem. In this case the factorization property, valid only in cartesian coordinates, is an 
exact numerical property and provides a strong stability test [|T|. 

The numerical methods used are based on the Hermite spline expansion, orthogonal 
collocation and iterative procedures for solving the linear system. An important step in 
its solution is done by means of the tensor trick M . 



3. RESULTS 

3.1. S-wave approximation 

We first discuss the results obtained in the S-wave approximation. By this we mean i) 
an interaction which acts only on S-wave and ii) that all the angular momenta variables 
in expansion (|3|) are zero. The interaction considered is MT I-III with parameters given 
in [[3], slightly different from those of its original version Q. The scattering length and 
effective range in different (T,S) channels are given in Table [l] and compared with the 
existing calculations. Our scattering length values are in agreement with || except in the 
S=l T=0 channel. For T=l they are also close to those obtained in || although with the 
original potential parameters. 

3.2. n-t scattering with realistic interactions 

A step behind the S-wave approximation has been achieved by considering the n-t 
scattering with realistic NN potentials. This reaction is a pure T=l channel, free from 
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Table 1 

Low energy N-NNN parameters Table 2 



T 


s 


a 


r 


a| 


4 He binding energy and 


r.m.s. 


radius 








14.75 


6.75 


*14.7 


'He i 





4 He* 





1 


3.25 


1.82 


*2.8 


B 30.30 30.31 


30.29 


0.27 


1 





4.13 


2.01 


4.0 


r.m.s. 1.44 




4.99 


1 


1 


3.73 


1.87 


3.6 









Table 3 

Experimental n-t scattering lengths 
fin fm.). 



a i 



Table 4 

Calculated n-t scattering lengths, zero-energy 
cross section and 4 He binding energy (B^) and 
r.m.s. radius (units are fm and MeV) 



3.91 ±0.12 


3.6 ±0.1 [12] 






a + 


ai+ 


a(0) 


B, 


r.m.s. 


4.98 ±0.29 


3.13 ±0.11 [ 






AV14 


4.31 


3.79 


193.7 


23.34 


1.56 


2.10 ±0.31 


4.05 ±0.09 
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Nijm II 


4.31 


3.76 


191.6 


23.39 


1.54 


4.453 ±0.10 


3.325 ±0.016 [14j 


MT I-III 


4.10 


3.63 


177.0 














exp 






170.2 ±9.7 







the difficulties related to the Coulomb interaction and for which accurate low energy 
scattering data exist pT| . The extrapolated zero energy cross section is cr(0) = 170.2 ±9.7 



fm 2 but the spin-dependent scattering length values are still controverted (see Table ^). 

We have calculated the scattering length of the singlet + and triplet 1 + states with the 
AV14 || and Nijmegen II |TI| potentials. The interaction is limited to the 1 S'o, 3 Si, 3 D\ 
waves and the FY expansion includes all the components with l z — 0, 1 and l y = 0, 1,2. 
The results are displayed in Table |4] together with the 4 He binding energy (B4) and r.m.s. 



radius. They are close to those presented in this Conference by the Pisa group |15| and 
our i? 4 value for AV14 shows a relative agreement of 10 -4 with |7j]. 

The MT I-III results presented in the previous section give a sligthly overestimated value 
for the zero energy cross section cr(0) = 184.6 fm 2 . However the inclusion of higher partial 
waves in the FY expansion reduces this value to cr(0) = 177.0 fm 2 in closer agreement 
with experiment. The corresponding scattering lengths are then a® = 4.10, a\ = 3.63 fm. 

The conclusion is that the considered realistic interactions on one hand give results 
very close to each other and on another hand fail in describing the zero energy n-t scat- 
tering as they failed in the 3N case. These problems are strongly correlated and can 
be solved by including a three- nucleon interaction (TNI). We have considered a simple 
phenomenological TNI model 

2w p-W 



W{p) = W r - W a - p=Jx 2 + y 2 

P P v 

and get the overall agreement B% = 8.48 MeV, B4 = 28.56 MeV, ao = 4.0 fm, a\ = 3.6 fm 
and o-(0) = 172.4 fm 2 with the parameters W r = 500 MeV, W a = 174 MeV and fi = 2.0 
fm . This set of values is far from being unique and the only serious constraint is the 
existence of a repulsive short range part. 
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The calculations have been completed up to the first 3-body break-up threshold by 
including the first negative parity states J n = 0~, 1~,2~ corresponding to a n-t relative 
P-waves. This region (T[ ab rs 3 MeV) is of particular interest due to the observed structure 
which has been interpreted as a series of P-waves resonances [|T(J . 

Let us first consider the results obtained with MT I- 1 1 1 , displayed in Figure [l]. This 
model conserves separately L and S and consequently the J n = 0~, 1 _ ,2~ states coming 
from S=l are degenerated. The corresponding cross sections differ only by statistical 
factors. The remaining J 71 = state comes from an (S = 0, L = 1) coupling. Despite the 
simplicity of this model the agreement with experimental data, specially in the resonance 
region, is very good. Only the zero energy region is a little overestimated. The n-t P- 
waves resonances are thus generated by S-wave NN interaction alone. The effective 1+3 
potential is essentially due to the exchange mechanism between the four nucleons. We 
notice however that a first attempt in describing this cross section with l z = 1 in the K 
components but keeping zero all the remaining angular momenta in ([3D failed. 




0.1 1 0.1 1 

E lab. (MeV) E lab. (MeV) 

Figure 1. n-t cross section (MT I-III). Figure 2. n-t cross section (AV14). 

The results obtained with the AV14 are plotted in Figure |2|. They are not able to repro- 
duce the experimental data neither in the S-wave nor in the resonance region. Correcting 
the zero energy with our TNI model results in lowering the value of the cross sections 
and the disagreement becomes worse. When compared to the MT I-III results one sees 
that, due to non central terms in the NN potential, the different J n contributions are 
shifted. As a consequence the coherent enhancement is weakened and the resonance pick 
broadened. On another hand the 2~ partial cross section, which should provide the domi- 
nant contribution, is noticeably smaller. Trying to account for this disagreement we have 
included the 1 Pq, 1 Pi, 3 P\ waves in Vnn, closing the J = 1 shell. The results, displayed 
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Table 5 

P-wave n-t phase shifts at E=3 MeV (degrees). 

Vjvjv partial waves 0~ 1~ 1~ 2~ 

AV14 1 S , 3 S U 3 D 1 21.2 26.6 47.8 35.4 

1 S , 3 S 1 , 3 D 1 + 1 P , 1 P u 3 Pi 24.1 20.0 37.5 32.7 

MT I-III ^o, 3 ^! 44.1 30.3 44.1 44.1 



in Table |5], show that the corrections go in the bad direction. On the contrary if the 
3 Si — 3 Di coupling is decreased and Vas 1 enhanced, keeping constant the triton binding 
energy, the experimental values are reproduced. 

The results obtained with MT I-III show the difficulty to understand these resonances 
in terms of the NN interaction alone. Several reasons can be put forward to explain the 
apparent difficulty for realistic interactions in describing the n-t scattering before having 
to conclude about its definite failure. Complementary calculations would be suitable with 
increasing number of FY components, by adding higher partial waves in the potential, by 
including realistic TNI, etc. They are however in the limit of the capacities of the best 
existing computers. 

On another hand it is intriguing that very simple interactions work so well. Is it only 
a chance that the conjugated actions of very sophisticated forces V 2 n + V^n + • • • result 
in a very simple two-body effective potential ?. Or should that rather be an indication 
that simple effective potentials must be introduced from the very beginning in these low 
energy few-nucleons physics ?. 
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